Strongly anisotropic non-equilibrium phase transition in Ising models with friction 
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The non-equilibrium phase transition in driven two-dimensional Ising models with two different 
geometries is investigated using Monte Carlo methods as well as analytical calculations. The models 
show dissipation through fluctuation induced friction near the critical point. We first consider high 
driving velocities and demonstrate that both systems are in the same universality class and undergo 
a strongly anisotropic non-equilibrium phase transition, with anisotropy exponent 8 = 3. Within 
a field theoretical ansatz the simulation results are confirmed. The crossover from Ising to mean 
field behavior in dependency of system size and driving velocity is analyzed using crossover scaling. 
It turns out that for all finite velocities the phase transition becomes strongly anisotropic in the 
thermodynamic limit. 
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I. INTRODUCTION 

The interest in magnetic contributions to friction due 
to spin correlations has strongly increased in recent years. 
One interesting aspect is the energy dissipation due to the 
formation of spin waves in a two-dimensional Heisenberg 
model induced by a moving magnetic tip [IH3 > which can 
be of Stokes or Coulomb type depending on the intrinsic 
relaxation time scales [I]. On the other hand, magnetic 
friction occurs also in bulk systems moving relative to 
each other. Kadau et al. [S] used a two-dimensional Ising 
model, cut into two halves parallel to one axis and moved 
along this cut with the velocity v, to explore surface fric- 
tion. The motion drives the system out of equilibrium 
into a steady state, leading to a permanent energy flux 
from the surface to the heat bath. This model exhibits 
a non-equilibrium phase transition, which has been in- 
vestigated in several different geometries [B] by means of 
analytical treatment as well as Monte Carlo (MC) sim- 
ulations. The critical temperature T c of the considered 
models depends on the velocity v and has been calculated 
exactly for various geometries in the limit v — > oo. In this 
limit the class of models show mean field-like critical be- 
havior. Subsequent investigations have been done in a 
variety of context, in particular for driven Potts models 
[7] and for rotating Ising chains of finite length [5] . 

The nature of non-equilibrium phase transitions is still 
a field of large interest, and simple models helping to 
explore this field are seldom. A very famous example is 
the driven lattice gas (DLG) f9iUl|. exhibiting a strongly 
anisotropic phase transition. Despite a lot of similari- 
ties between the driven lattice gas and the Ising model 
with friction, there is an important difference: The order 
parameter is conserved in the former, while it is non- 
conserved in the latter model. A further class of models 
characterized by non-equilibrium phase transitions are 
sheared systems |I2lTH] . experimentally accessible within 
the framework of binary liquid mixtures. 

Like the driven lattice gas, the systems investigated in 
the following exhibit a strongly anisotropic phase tran- 
sition, which is investigated by means of Monte Carlo 



simulations as well as a field theoretical ansatzln addi- 
tion, the case of finite velocities v is analyzed by means of 
crossover scaling, where a broad range of velocities and 
system sizes are analysed. We show that for all v > the 
considered models end up in the mean field class with 
strongly anisotropic correlations as soon as the system 
size exceeds a velocity dependent crossover length L x (v). 

While a crossover behavior from Ising to mean field 
class occurs in various thermodynamic systems such as 
ionic fluid |15III6| and spin systems with long-range inter- 
actions [T7j, to our knowledge such a crossover including 
a change from isotropic to strongly anisotropic behavior 
has not been investigated in detail until now. The pa- 
per is organized as follows: After introducing the model 
and geometries, we determine the anisotropy exponent 
for v = oo using MC simulations as well as a field theoret- 
ical model. Then we turn to finite velocities and present 
the crossover scaling analysis. Finally we discuss our re- 
sults. 



II. MODELS 

The systems considered in this work are denoted 2d 
and 1+Id and are shown in Fig. [I] (for a classification 
see [1]). The 2d system is a two-dimensional two-layer 
Ising model with Lu x L± x 2 lattice sites, where the two 
layers are moved relative to each other along the paral- 
lel direction. Each lattice site carries one spin variable 
dij^ = ±1, and only nearest-neighbor interactions are 
taken into account. Periodic boundary conditions are ap- 
plied in both planar directions, i.e., <Ji.j t k = <^i+L^,j,k = 
&i,j+L ± ,k- In order to simulate a finite velocity v using 
Monte Carlo simulations the upper sub-system is moved 
v times by one lattice constant during each random se- 
quential Monte Carlo sweep (MCS). Since one MCS cor- 
responds to the typical time to « 10~ 8 s a spin needs to 
relax into the direction of its local Weiss field, and as the 
lattice constant is of the order ao ~ 10~ 10 m, the velocity 
v is given in natural units ao/tg ss lcm/s. 

Instead of moving the two layers against each other, 
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Figure 1. (Color online) The systems considered in this work. 
The gray regions represent the magnetic systems, while the 
green (dark) regions are the moving boundaries. The arrows 
indicate the motion of the subsystems. 
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Figure 2. (Color online) Cross sections of the 1+ld (a) and 
the 2d model (c), and slight modifications of both models 
((b) and (d)). The grey circles represent spin chains and 
the connecting lines substitute for the coupling, where green 
wiggled lines stand for moving and black lines for stationary 
couplings. Black crosses and dots indicate a motion into and 
out of the plane, respectively. 



we reorder the couplings between the subsystems with 
time to simplify the implementation [5J. Introducing a 
time-dependent displacement 

A(t) = vt, (1) 

which is increased by one after each 2LnL±/v random 
sequential spin flip attempts, the Hamiltonian can be ex- 
pressed as 

i=i;=i*=o 

L\\ L ± 

%=\ j=i 

with the reduced nearest neighbor coupling K — ft J , the 
reduced boundary coupling — /3Jb, and /3 — l/k^T. 
In the following we assume J = Jb — 1. 

The critical temperature T c (v) of the regarded systems 
increases with v and saturates for high velocities. In the 
limit v — > oo an analytical calculation of the critical tem- 
perature for the 2d geometry yield 

T c 2d (oo) = 4.058782423... (3) 

for J = Jb = 1 [6 . The basic idea of the analytic solution 
provides the approach for the implementation of infinite 
velocity, which works as follows: the interaction partner 
for a spin in the lower layer is chosen randomly from the 
same row in the upper layer. Thus we can use Eq. |2| 
with a random value 1 < A(t) < Ln. 

The 1+ld system consists of a two-dimensional Ising 
model, where all rows are moved relative to each other. 
The displacement A(t) = vt as well as the coupling K± 
is equal for all adjacent rows, leading to the Hamiltonian 

L|, L ± 

pH{t) = - ^2^2Kna i ja i+hj +Kx(Ti, j (7 i+A(t)>j+1 . (4) 

<=1 3=1 

Again, periodic boundary conditions are applied in 
both directions, where discontinuities in _L direction are 



avoided through the homogeneous displacement A(t) |6_. 
The analytical treatment at v — > oo gave the critical tem- 
perature 

T c 1+ld (oo) = l/log^A/S + vTr) = 3.46591... (5) 

for J|| = Jj_ = 1 in this case [BJ- Within the scope of 
the 1+ld model the velocity v corresponds to a shear 
rate, which is ofter denoted as 7 [181 IH)] . However, we 
will use the term velocity for both driving mechanisms 
throughout this work. 

In the following we argue that both systems show the 
same underlying critical behavior. In order to emphasize 
the similarity, Fig. [2] illustrates slight variations of both 
models. First of all we start with the 1+ld model (a) 
and change every second bond perpendicular to the mo- 
tion into a stationary bond. Additionally, we perform a 
transformation that changes the homogeneous shear A(t) 
into an alternating shift ±A(i) of the double chains and 
reverses (i — > —i) every second double chain, leading to 
the configuration in Fig. (2Jd. These modifications do not 
change the critical behavior of the 1 + ld system, since 
still one-dimensional chains (now consisting of two rows) 
are moved relative to each other. On the other hand, 
the cross section of the 2d model can be visualized in a 
slightly different way (see Fig. I2H) without altering the 
corresponding Hamiltonian, Eq.T2j). Since the next near- 
est double chains in (b) are not moving relative to each 
other, the only difference between (b) and (d) are the 
third nearest neighbor bonds in (d), which are irrelevant 
at the critical point where long range correlations dom- 
inate. Hence we conclude that both systems belong to 
the same universality class. 

Finally we mention that we must use the multiplicative 
rate 

p nip (AE) = e-^ AE - AE ^, (6) 

with A-E m i n = min({AS}) to reproduce the critical tem- 
peratures, Eqs.([3] [5]), in simulations, for a discussion see 
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Figure 3. Snapshots of one layer of the 2d model with Lu = L± — 512 and J« = J± = 1 at temperature T = 3.5, which fulfills 
T c 2d (0) < T < T c 2d (oo). We start with an equilibrium system at v = (left), set v — oo, and show the evolution at t — 42 MCS 
(center) and t = 360 MCS (right). 



III. RESULTS 

In order to illustrate symptomatic features of both sys- 
tems, Fig. [3] shows a sequence of spin configurations of 
one layer of the 2d system (note that the same char- 
acteristics are observed in the 1+ld system). On the 
left hand side an equilibrated system at T = 3.5 well 
above the critical temperature of the non-moving sys- 
tem, T c 2d (0) = 3.20755(5) gU], is presented. Shortly after 
starting the motion stripe-like domains arise, spanning 
the whole system parallel to the motion. The stripes are 
rather stable, but are nonetheless transient, since they 
grow in time until the system ends up in a homogeneously 
magnetized state. The evolution in Fig. [3] is an example 
for a velocity driven phase transition already described 
in [6] , which is triggered by the onset of the motion 
and the associated increase of the critical temperature. 
The circumstances are comparable to a quench, which is 
characterized by a temperature decrease below T c . After 
a quench a coarsening of domains is observed, whereas 
the growth of the domains can be described by a power 
law (e.g. |211 |2"2"] ). Domain growth in systems exhibiting 
a strongly anisotropic phase transition, e.g., the DLG 
model, is also a well investigated subject |10[ |23"1 124] , 
The corresponding time evolution of spin configurations 
are similar to those shown in Fig. [3] leading to the as- 
sumption that the 2d and the 1+ld geometries are also 
characterized by strongly anisotropic correlations, which 
is shown in the following section. 



dependent critical exponents at the critical point |25| , 

Ut) ~°lp*-^ s (7) 

with direction fi — {_L, ||} and reduced critical temper- 
ature t = T/T c — 1. Defining the anisotropy exponent 

[HH2E] 



we find 

£n(*)/&(t)~£y& (9) 

independent of t. Isotropic scaling takes place for 
9=1 and strongly anisotropic scaling is implied by 
5^1. Several models with strongly anisotropic be- 
havior where studied in the past. Examples are Lif- 
shitz points as present in the anisotropic next near- 
est neighbor Ising (ANNNI) model [2HJ [30], the non- 
equilibrium phase transition in the DLG |10| . the two- 
dimensional dipolar in-plane Ising-model [28J. Further- 
more, strongly anisotropic behavior usually occurs in dy- 
namical systems, where the parallel direction can be iden- 
tified with time and the perpendicular direction(s) with 
space |27l |3"T] . In the latter case the anisotropy exponent 
9 corresponds to the dynamical exponent z. 

The knowledge of the anisotropy exponent is essen- 
tial and necessary for appropriate simulations of strongly 
anisotropic systems. To avoid complicated shape effects 
it is required to keep the generalized aspect ratio |2"51 - f2"8] 



A. Determination of 8 in the limit v — > 00 n = — qq~\ 

A strongly anisotropic phase transition is characterized fixed, which requires the knowledge of 9. We will show in 
by a correlation length £ M which diverges with direction the following that the limit p — ¥ simplifies the analysis 
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Figure 4. (Color online) Rescaled correlation function G± (Lm , oo; r±) at criticality for both models for varying system extensions, 
L|| = {40,80,160,320,625,1250,2500,5000,10000} (1+ld) and Ly = {16,32,64,128,256,512} (2d), respectively. The insets 
show £x(L|| ) whereby we yield £± by fitting an exponential function to the long-range part of G±(L\\ , oo; r±). The solid line is 
a power law with exponent _1 = 1/3 as predicted by the field theoretical analysis (see text). 



for infinite velocity v and turns out to be essential at 
finite v. 

We first discuss the case v — ¥ oo and always assume 
criticality, t = 0. In order to determine the anisotropy 
exponent 6 we calculate the perpendicular correlation 
function Gi_(L\\ , Lj_; rj_) = (ai_jCTi.j +r± ) between spins 
at distance r± in cylinder geometry L± — > oo (leading to 
p — > 0), and thereby gain the correlation length 
through 



Gl(£||,oo;7_l) ~ d±(L\ l )e- rx ^ L i\ 



(11) 



where the prefactor G±(L^) is shown to be proportional 
to L|| 2 ^ 3 in Appendix |a Approaching the critical point 
within the given geometry, the correlation length £y (t) is 
limited by Ly , and using Eq. (|9j) this leads to the relation 



U(h)~A ± Li 



(12) 



with non-universal amplitude A± [28 , 32J . Measuring the 
correlation length £j_ in dependency of the parallel exten- 
sion iy allows us to determine the anisotropy exponent 
9. 

In the simulations, the limit L± — >• oo is imple- 
mented by the condition L±/£± > 10. This is suffi- 
cient to keep the sytematic errors in G± smaller than 
the statistical error e = 1CP 3 adequate to calculate 
From e we can determine the required system sizes via 
L±/£± = — 21og[e/G_i_(L||)], where the factor 2 accounts 
for the periodic boundary conditions. As G± « 0.1 for 
L|| = 40 and G x ~ 0.02 for L\\ = 10 4 for the 1+ld model 
(see Fig. gleft)) we yield « 10 for L\\ = 40 and 

L±/£± ~ 0.7 for L|| = 10 4 , meaning that for large sys- 
tems a much smaller value of L±/£± would be sufficient. 

Fig. [3] displays the correlation functions for both mod- 
els. For the 1+ld case these correlations are purely ex- 
ponential also at short distances, since the coupling in 



+ direction is mediated through fluctuating fields [3], 
leading to dimensional reduction to an effectively one- 
dimensional system. The resulting correlation length £j_ 
is shown in the inset of Fig. |4^1eft) . The growth of £j_ {L\\ ) 
follows a power law with exponent = 1/3 and with 
prefactor 

A 1 + Xd = lim 2+ 1/3 £i +ld (L||) =0.68(2), (13) 



indicated as a black line. 

In the case of the 2d model (right figure in Fig|4j 
we find two regions with different characteristics. The 
short-distance correlations are affected by the + nearest- 
neighbor interactions within the planes which are not 
present in the 1+ld model. These correlations decay 
with a correlation length of the order ^ q [T c 2d (oo)] « 1. 
For large distances the correlations crossover to an expo- 
nential behavior. The exponential correlations are prop- 
agated by the fluctuations of stripe- like domains. The 
analysis yields 



Af= lim 



" 1/3 £i d (£||) = 0.94(3) 



(14) 



in this case. 

From the anisotropy exponent 6 = 3 we can derive 
the correlation length exponents v\\ =3/2 and v±_ = 1/2 
using the generalized hyper-scaling relation 



2 - a = 2/3 + 7 = v« + (d - l)u ± , 



(15) 



with d = 2 and mean field exponents a = 0, /3 = 1/2, 
and 7 = 1, whose validity has been demonstated in [B] 
by a mapping onto a mean field equilibrium model. 

The calculatation of 9 in the limit v — > oo is 
done within a one-dimensional Ginzburg-Landau- Wilson 
(GLW) field theory [33]. For v — >• oo it was shown in 
Ref. [B] that the 1+ld model can be mapped onto an 
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equilibrium system consisting of one-dimensional chains 
that only couple via fluctuating magnetic fields. Due to 
the stripe geometry with short length Ln and the periodic 
boundary conditions in || direction the magnetization is 
homogeneous in || direction, and parallel correlations are 
irrelevant. Hence we can use the zero mode approxima- 
tion in this direction. However, it is necessary to include 
a term representing the interaction between adjacent spin 
chains. This can be expressed by the square of the spa- 
tial derivative of the magnetization in the direction _L to 
the motion. Hence the minimal GLW model to describe 
this strongly anisotropic mean field system is given by 



PH = Ln 



,lr I ^m(x) 2 + ^m'(x) 2 + ^m{x) 4 ) (16) 



with phenomenological parameters t and u, where m(x) 
represents the magnetization of the spin chain at _L coor- 
dinate x. Eq. (16 I corresponds to the Hamiltonian used 



for the description of a cylinder-like spin system, which 
is infinite along one dimension, and finite and periodic in 
d— 1 dimensions [33]. The partition function of Eq. (16 1 



can be mapped onto a one-dimensional Schrodinger equa- 
tion in a quartic anharmonic oscillator potential using a 
rescalation, which yields the critical exponents m =3/2 
and 9 — 3. The detailed derivation is given in Ap- 
pendix [A] 



B. Crossover scaling at finite velocities 



Fig. [5k shows the unsealed data, which gives evidence 
that the correlation length of systems moved at high 
velocities v are well described by the exponent 9 = 3 
(dotted line), whereas for low velocities v < 2~ 12 ef- 
fectively the Ising exponent 9 = 1 (dashed line) holds 
for the simulated system sizes Ln. The curvature of 
the data of intermediate velocities suggest the crossover. 
As a data collapse on the analytical known [33] relation 
£_l(0,0,Z||)/L|| ~ ^4 oq = 4/tt (dashed line in Fig. ^ has 
to be obtained in the limit v — > 0, both axes must be 
rescaled by the same factor L x (v). This crossover length 
can be determined by applying the following method: We 
start with plotting the correlation length in the mean 
field limit £j_(t c (oo), oo, Ly). Then we subsequently add 
the data for smaller v by rescaling £j_ and Ly with L x 1 ( v ) > 
which shifts the points parallel to the dashed line, until 
a data collapse is obtained (see Fig. [5b) . This procedure 
works quite accurate for velocities v~> 2 -3 , only at very 
small v < 2~ 12 the errors in L x (v) grow due to the fact 
that we just shift the data along the dashed line. The re- 
sulting crossover length L x (v) is pictured as black dots in 
the inset of Fig. [5^,. The behavior of L x (v) is analogous 
to the velocity dependency of other quantities like the 
critical temperature or the energy dissipation, which are 
characterized by a power law for tiCl and a saturation 
for »>1. 

We conclude that for all finite velocities v > the crit- 
ical behavior changes from Ising type to mean field type 
at a velocity dependent crossover length L x (v) approxi- 
mately given by 



We now turn to finite velocities. The following analysis 
is exemplarily done for the 1+ld model, but as stated 
above, both models belong to the same universality class 
and similar results are expected for the 2d model. As 
we expect a crossover from an isotropic Ising model with 
9 = 1 to a strongly anisotropic system with 9 = 3, we 
must be careful with the system geometry: We cannot 



use a fixed finite generalized aspect ratio p, Eq. (10 1, in 



the simulations, as 9 is not constant. The only possible 
choice is p — > (or p — ¥ oo), where the ^-dependency 
drops out. 

We consider the correlation length £±(t c {v), v, L») at 
reduced critical temperature 



t c (v) 



Tc(v) 
Tc(0) 



1. 



(17) 



where T c (0) = 2/log(V2+ 1). t c (v) is calculated via 
a finite-size scaling analysis of the perpendicular corre- 
lation length (not shown). As this procedure becomes 
inaccurate for small velocities v < 2~ 8 , we calculate the 
critical temperature according to 



t c (v) v ^°cv<t 



(18) 



with c = 0.29(1) in these cases, where we assume <fi = 1/2 
in agreement with the literature [TU1 IT51 HU] , The results 
are shown in the inset of Fig. [5k. 



L x (v) 



3/2 



(19) 



(solid red curve in the inset of Fig. [5k), where the velocity 



is measured in units 10 m/s and the size in 10 m. 
The velocity independent prefactor was added to shift the 
crossover point, i.e., the intersection of the asymptotes, 
to z = 1. The saturation of L x at v x = 18(2) results 
from the lattice cut-off, as L x (v x ) rj 1. The inset in 
Fig. [SJd shows the effective exponent 9 e g, obtained from 
the logarithmic derivative 



y cff 



d log L II ' 



(20) 



whose value changes from 9 c g = 1 (Ising, isotropic) to 
#eff = 3 (MF, strongly anisotropic). Note that we veri- 
fied the mean field exponents for v > 1/8 with finite-size 
scaling methods and also found good agreement of the 
scaling function with the universal finite-size scaling func- 
tion [35] (not shown). In order to illustrate the change of 
the critical behavior, Fig. [6] shows typical critical spin 
configurations for different values of the crossover scaling 
variable z = L\\/ L x {v). 

We are now able to compare our results with the lit- 
erature. If the crossover scaling variable z<l Ising-like 
behavior occurs, whereas for z> 1 mean field exponents 
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Figure 5. (Color online) Velocity dependent crossover behavior in the 1+ld case. Both pictures show log-log plots of the 
correlation length £x(t c (v), v, L«) as function of the system size Lii at reduced critical temperature t c (v) for a broad range of 
different velocities v. The dashed line is the analytically known Ising limit £x(0, 0, L^)/L^ ~ 4/-7T valid for v — > [34 , while 
the dotted line has slope (9 _1 = 1/3. The left figure shows the unsealed data and the inset displays the rescaling factor L x (v) 
for different velocities v (black dots, see text) and a function approximating the data given in Eq. (19 1 (red solid line), as well 
as the reduced critical temperature t c (v) (blue squares) together with its asymptotes, Eqs. ( |5|18| |. The right figure displays the 
same data rescaled with the crossover length L x (v), leading to an excellent data collapse. The inset shows the crossover of the 
effective anisotropy exponent 8 c g from 9 c g — 1 (Ising, dashed line) to o g = 3 (MF, dotted line). 



and strongly anisotropic correlations are expected. In ex- 
periments [13 , even slow shear rates of the order of 10 -4 
(in natural units ijj" , where now to is the time scale of 
the fluid dynamics), lead to a crossover length L x < 100 
and, as the typical system size is large wrt. the atomic 
distances, give z» 1, indicating that experimental data 
are always obtained in the mean field limit. 

In relation to the results of Winter et al. [19] we find 
that the correlation length exponent has been measured 
in the regime 29 < z < 239, leading to the anisotropy ex- 
ponent « 3 in agreement with our results. In Ref. [18 
the correlation length exponents have also been deter- 
mined in the mean field limit. Looking at the lowest 
velocity v = 1/32 we find 53 < z < 1066, where a sur- 
prisingly small anisotropy exponent 9 » 0.73 has been 
estimated. The highest velocity v = 50 leads to 9 « 1.2 
and 1100 < z < 22000. These discrepancies might be 
attributed to the fact that an integral quantity, the or- 
der parameter, has been measured, as well as to strong 
surface effects induced by the open boundary conditions 
used in the _L direction. 



IV. CONCLUSION 

In this work we investigated two recently proposed 
driven Ising models with friction due to magnetic in- 



teractions, namely the 1+ld and 2d model, using MC 
simulations as well as analytical methods. At first we 
focused on the strongly anisotropic critical behavior and 
calculated the anisotropy exponent 9 in the limit of high 
driving velocity v — ¥ oo. Therefore the perpendicular 
correlation function of a cylinder-like geometry was cal- 
culated at criticality for different system sizes. Evaluat- 
ing the connection between system size and correlation 
length, Eq. (12 1, we were able to find the critical expo- 
nents 9 = 3 as well as v\\ = 3/2 and v± = 1/2. The 
analytic deviation of these exponents within the frame- 
work of a Ginzburg-Landau- Wilson Hamiltonian led to 
the same values. Comparing the results to the driven 
lattice gas [HI [TU] we note that it also shows a strongly 
anisotropic phase transition at a critical temperature 
which grows with the velocity. Remarkably this phase 
transition is characterized by the same critical exponents 
at large fields. 

Finally we focused on the critical behavior for finite 
velocities v and performed extensive MC simulations in 
order to calculate the crossover scaling function describ- 
ing the crossover from the Ising universality class at v = 
to the non-equilibrium critical behavior at v — > oo. The 
analysis has exemplarily been done for the 1+ld model, 
but as shown, both models belong to the same univer- 
sality class and similar results are expected for the 2d 
model. In the analysis an additional complexity arised 
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Figure 6. Typical spin configurations of the critical 1+fd system for L« = 64 and different velocities v = 2 , . . . , f . z = 
L\\/L x {v) denotes the crossover scaling variable (see text). The critical domains are isotropic and Ising-like for z < 1 and 
become anisotropic for z > 1. 



due to the strongly anisotropic characteristics of the cor- 
relations. Therefore we calculated the correlation length 
in a cylindrical system, circumventing intricate shape ef- 
fects. We were able to identify a crossover length L x (v) 
using a simple method based on the rescaling of data for 
each velocity such that a data collapse occurs. This pro- 
cedure leads to an excellent data collapse of all simulation 
results for different velocities v and system sizes Ly. 

It turns out that for all finite velocities v > the mod- 
els undergo a crossover, at crossover length L x (v), from 
an quasi-equilibrium isotropic Ising-like phase transition 
to a non-equilibrium mean-field behavior with strongly 
anisotropic correlations. 
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Appendix A: Scaling exponents of the GLW model 

The following calculation is similar to [33 . Discretizing 
the integral 

m = L,| £dx ( *m(z) 2 + ^m'(x) 2 + ^rn{x)^j 

(Al) 



with step size Sx, NSx = L±, m, = m{i5x) and 5nii = 
m i+ i - mi gives 

Pn = L ]l f:sJ t -^ + 1 - 5 ^ + ^ m t). (A2) 
i=i ^ ' 

In order to evaluate the partition function 

/oo 
V[m(x)]e-^, (A3) 
-oo 

we use abbreviations in analogy to transfer matrices, 

reft 2 . u 4 \ l~T L,,Sm 2 

T(m,m+)=e- L "MF»+4l™) ^^"4^, (A4 ) 

s v / s ✓ 

V(m) U(5rn) 

with m + = m + Sm to get 

/oo /*oo *>co 

dmi / dm2T(mi,m 2 ) / dm 3 T(m 2 , m^) x 
-oo J —oo J —oo 

/oo 
dmjvT(mAr_i, m7v)T(mAr, mi) (A5) 
-oo 

for the assumed periodic boundary conditions. 

Let x + = x + Sx and ip(m + ) be the result of the in- 
tegrations for the interval ]a; + ,LjJ. Since T(m, m + ) is 
near-diagonal for Lu — > oo, we can write -0(m + ) as 

X^(m + ) fa ip(m) + ip'(m)Sm + \^)" \m)8m 2 , (A6) 

where A denotes the growth factor of the integrations 
corresponding to the leading eigenvalue of the transfer 
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matrix T(m,m + ). The integral over m + in the partition 
function becomes 

/>oo 

i/j(m) = / dm + V(m)U(m + — m)i/)(m + ) 



V(m) 

V(m) 
A 



>o L||(5m 2 



^M + ^''(m)), (A7) 



and yields the solution of the integrations for the interval 
[x + , LjJ. Hence we get a differential equation for ip(m), 

V(m) (ip(m) + ^-/(m)) = A^(m). (A8) 



We now substitute 

ipirn) — > 'J'(m) 



m -> mu- 1/6 L~ 1/3 
A -> 1 - A<fccit 1/3 Ljj~ 1/3 



(A9a) 
(A9b) 
(A9c) 
(A9d) 



and expand to lowest order around L\\ — oo to yield the 
Schrodinger equation in a quartic potential, 



-^ + |m 2 +^m 4 -A)*, / ;, 1 =0. ,A1D] 



valid in the scaling limit — ^ oo, t — > with x = 
t{L\\lu) x l v w kept fixed. 

The correlation length £_l(£m) is determined from the 
lowest eigenvalues Aq,i of this equation, as 



£_l = 5x[ log 



1 



A, 



Ai / Ai — Ao V u 



1/3 



(All) 



From the substitution, Eqs. ( A9 ), we directly read off the 



exponents vu = 3/2, and 9 = 3. 

The correlation function amplitude G±(Lu) from 



Eq. (Ill is proportional to m 2 and thus scales as L, 



-2/3 



as canoe seen from Eq. ( A9b ) 
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